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Purpose: Beam range control is the essence of radiotherapy with heavy charged 
particles. In conventional broad-beam delivery, fine range adjustment is achieved 
by insertion of range shifting and compensating materials. Ideally, such material 
should be water equivalent as well as that for dosimetry. In this study, we evaluated 
dosimetric water equivalency of four common plastics, HDPE, PMMA, PET, and 
POM, by uniformity of effective densities for carbon-ion-beam interactions. 
Methods: Using the Bethe formula for stopping, the Gottschalk formula for multiple 
scattering, and the Sihver formula for nuclear interactions, we calculated the effective 
densities of the plastics for these interactions. We tested HDPE, PMMA, and POM 
in carbon-ion-beam experiment and measured attenuations of carbon ions, which 
were compared with empirical linear-attenuation-model calculations. 
Results: The theoretical calculations resulted in reduced multiple scattering and 
increased nuclear interactions for HDPE compared to water, which was verified in 
the experiment as increased attenuation of carbon ion by 8.9% for 20-cm range shift. 
The increase of attenuation was 2.5% for PMMA and 0.0% for POM for similar 
range shift. PET was theoretically estimated as in between PMMA and POM. The 
agreement between the measurement and the calculation was typically about 1%. 
Conclusions: For range control, POM was the best and virtually indistinguishable 
from water. PET was better than PMMA, which was also reasonably water equiva- 
lent. Large range shift with HDPE should be avoided for increased attenuation. The 
simplistic nuclear-interaction model agreed with measurement at the level of 1% in 
attenuation of carbon ions. 
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I. INTRODUCTION 



The essence of radiotherapy with heavy charged particles is intrinsic capabihty of three- 
dimensional dose formation with a Bragg peak at the beam range, which can be precisely 
controlled by choosing the beam energy incident on a patient. In broad-beam- delivery 
35 techniques,-*^ a range shifter reduces excessive beam energy and a range compensator com- 
pensates the deficits of a patient body so that the beam range in the field will match the 
distal surface of the planning target volume for the best dose conformation.- 

A range shifter is typically a composite of plates of precisely controlled thickness and a 
range compensator is a physical filter custom-made for an individual field. Their material 

40 should be uniform, stable, machinable, inexpensive, and ideally equivalent to water that is 
the reference material for dosimetry. These requirements are also valid for solid phantom 
materials that are used when water-phantom dosimetry is not possible or convenient. To es- 
timate the effective thickness of a material by the resultant range shift in water, an effective 
density is defined for the material regarding as variable-density water (H2O). However, this 

45 approximation may lead to dosimetric errors because the effective density differs for multiple 
scattering and nuclear interactions, which also influence the dose distribution.^ Dosimetric 
water equivalency may be evaluated by uniformity of effective densities for relevant interac- 
tions, which is essential for materials for range shifting, range compensation, and phantom 
dosimetry. 

50 Some of common materials for such purposes are high-density polyethylene (HDPE),- 
polymethyl methacrylate (PMMA),- and synthetic resin of acrylonitrile, butadiene, and 
styrene (ABS).- In general, HDPE is inexpensive and close to water in density, PMMA 
is better at dimension stability and contains oxygen that dominates water, and ABS is 
available in various forms as chemical wood in industry. For proton and ion-beam dosimetry, 

55 PMMA has been commonly used as water-equivalent phantom material.-i^ In this study, we 
quantitatively evaluate dosimetric water equivalency of these plastics. 
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TABLE I. Material properties (elemental composition, oxygen weight fraction ws, typical density 
p, number of electrons per atomic mass unit {Z/Aj.), and mean excitation energy /) excerpted from 



Ref. M (also available at the Particle Data Group website http://pcig.lbl.gov/) except for an 
empirical HDPE density. 



Material Composition 


m 


P 


{Z/A) ^ 




HDPE (C2H4)„ 


0.000 


0.96 


0.57034 57.4 


PMMA (CsHgOs)™ 


0.320 


1.19 


0.53937 74.0 


PET (CioH804)„ 


0.333 


1.40 


0.52037 78.7 


POM (CHaO)^ 


0.533 


1.42 


0.53287 77.4 


Water H2O 


0.888 


1.00 


0.55509 79.7 



II. MATERIALS AND METHODS 



A. Sample plastics 

In addition to HDPE and PMMA, we theoretically evaluated two other oxygen-rich plas- 
60 tics, polyethylene terephthalate (PET) and polyoxymethylene (POM), among which we 
experimentally tested HDPE, PMMA, and POM with carbon-ion beams. We excluded ABS 
as it varies in composition of carbon, hydrogen, and nitrogen. Table [T] summarizes the rel- 
evant properties of these plastics and water, ^ in which the densities may vary with degree 
of polymerization. Incidentally, electron density rif. = (Z/Ar) p/u is the effective density for 
65 Compton scattering that dominates the attenuation of megavoltage photons. 



B. Effective densities for ion-beam Interactions 

1. Stopping-power ratio 

The effective density for range shifting is the stopping-power ratio of material to water. 
In the Bethe theory, it is defined as 

S _ {Z/Ar)p In gHMc! _ 1^(1 _ _ g 

(Z/Ar)^ p^ln'-f^- ln(l ^ ^ 

J-w C C 
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where properties p, {Z/A^), and / for the plastics and those for water, p^,, {Z/Ar)w, and 
Iwi are given in Tabled me is the electron mass, and v and c are the particle speed and the 
light speed in vacuum. For materials with / ^ 1^, the f-dependence of S/S^ is small. For 
example, S/S^, varies by only 0.07% for full stopping in HDPE from f^/c^ = 0.5 or nucleon 
kinetic energy E/A = 386 MeV. As there is generally large uncertainty in the p and / values 
for compouds, the scattering-power ratio is usually measured by range shift s in water per 
material thickness t inserted upstream of water, 

S ^ {Z/Ar)p \n'-^ ^ s 

ignoring the amount of the air replaced by the material. 



2. Scattering-power ratio 

70 The effective density for multiple scattering is the scattering-power ratio of material to 
water. Using the scattering-power formula for heavy-charged particles by Gottschalk,^- it is 
defined as 

^ = y (29.73 -\nZi- In A,^) , (3) 

91.69p^ i ' '''' ^ ' 

where Zj, Ari, and Wi are the atomic number, the atomic weight, and the weight fraction of 
75 element i. 

For a certain range shift s = {S/Su])t, the increase of mean square angle of the primary 
particles in the material differs from that in water by a double-ratio factor, 

A^(^) = ^A^4^)- (4) 
The increase of mean square angle in water can be estimated by a semi-empirical formula,— 

1000 \^p/ ^ — s 

for particles of range Rq, charge ze, and mass {m/mp) in units of the proton mass. 
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3. Nuclear- cross- section ratio 



The effective density for nuclear interactions is the nuclear-cross-section ratio of material 
to water, 

w Pw 

where 

^^a/A^) = J2 ^^A. (7) 

is the nuclear cross section per atomic mass unit of the material and aAi is the collision cross 
section between a projectile nucleus of mass number A and a target nucleus of element i. 
Using the Sihver formula,— it is defined as 

r .1 . ' ' ^ '^'^ 



1 , -i 
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(8) 

2.247- 0.915(^-3 + 1) fori = l 
boi={ ^ ^ (9) 

1.581 - 0.876 (^A-^i + An ^) for i > 1, 

where tq = 1.36 fm is the effective nucleon radius and we applied the proton-nucleus formula 
to transparency parameter bo for collision on hydrogen (z = 1). The energy dependence is 
ignored as it is almost common to the relevant nuclei for collisions with E/A ^ 100 MeV in 
the Sihver model. 

Primary-particle fiuence $ decreases in matter due to nuclear interactions. With an insert 
of range shift s, the fiuence at depth d in water is factorized as 

^d) =a{s)%{s + d), (10) 

where ^o{d) is the original fiuence at depth d in water without an insert and a{s) = 
$(0)/$o('S) is the ratio of the fiuence attenuated by the insert to the original fiuence at 
depth s in water. The fractional fiuence loss per range shift, —d^/{^ds), differs between 
the material and water by a double-ratio factor. 



1 d^ cr/a^ f 1 d^\ 
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$ ds S/S^, \ ^ ds J , 
which is solved with boundary condition $ = $o for s = 0, resulting in 



(11) 



If we employ a simplistic linear-attenuation model,— 



(13) 



where relative attenuation power is the ratio of fiuence loss per depth in water to the 
final fiuence at the range, then the attenuation ratio for range shift s reduces to 



Matsufuji et al. measured attenuation of carbon ions in PMMA resulting in $(i?o)/'^'(0) ~ 

1 

0.47 for i?o = 13.5 cm as shown in Fig. 5 of Ref. il5|, which translates into ~ /3pmma = 
0.084 cm^^ if we assume PMMA to be sufficiently water equivalent. 

85 C. Experiment with carbon-ion beams 

The equivalency of material to water should be evaluated by equality of the resultant 
doses. As far as the primary particles are concerned, the Bethe theory for stopping and 
the Moliere theory for multiple scattering are accurate to the level of measurement limit if 
the p and / values are adjusted to specific material.—*^ On the contrary, there only exist 
90 experimental cross-section data and empirical models for nuclear interactions. Although the 
model chosen for this study has been generally tested,— specific validation may be needed 
for this particular application. 

An experiment was performed with carbon- ion beams of E/A = 290 and 430 MeV ex- 
tracted from Heavy Ion Medical Accelerator in Chiba (HIMAC) of National Institute of 
95 Radiological Sciences. The beams were laterally broadened by a wobbling/scattering sys- 
tem to form a uniform field of 10-cm diameter,— longitudinally moderated by a ripple filter 
for Gaussian range modulation of 1.8 mm (rms),^^ and delivered horizontally to a box-shaped 
water tank on a movable treatment couch. 

We used binary plates of HDPE, PMMA, and POM, which were stacked at the immediate 
100 upstream of the tank to form an insert of arbitrary thickness t in steps of about 1 cm. To 
eliminate the infiuence of the beam divergence by wobbling at about 10 m upstream, the 
treatment couch was moved downstream so that the Bragg peak would stay in the same 
place in the laboratory system. The amount of couch movement was the expected range 
shift s = {S/S^) t using our standard S/Su, values 1.02 for HDPE, 1.16 for PMMA, and 1.36 




(14) 
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20 40 

Horizontal position / cm 



FIG. 1. Side views of apparatus for carbon- ion beam experiment (a) without and (b) with a plastic 
insert of thickness t and couch movement by range shift s, where the central-axis (horizontal-dotted- 
line) doses at varied depth d were measured with PPIC. 



105 for POM. The central-axis doses at varied depth d in water after the tank wall of 19-mm 
PMMA were measured with a Markus-type plane-parallel ionization chamber (PPIC) with a 
protective cover of 0.87-mm PMMA (PTW Type 23343), as shown in Fig. [TJ The dosimetric 
precisions were as good as 0.1% for peak dose and 0.1 mm for depth in reproducibihty. 

We chose insert thicknesses that approximately corresponded to range shifts of 4 cm 
and 8 cm for 290 MeV and 10 cm and 20 cm for 430 MeV. For beam deliveries of equal 
monitor units, the central-axis doses were measured at minimum depth intervals of 0.1 mm 
around the moderated Bragg peak. For the reference depth of the moderated Bragg peak, 
we adopted depth d^Q, at which the dose is decreased to 80% of peak dose Dp. Range shift 
s was measured as the shift of dgo by insertion. Using the dose points in the depth region 
of dso + 1 cm ^ d < dso + 2 cm, we estimated the fragment dose -D/ by extrapolation of the 
linear- fit line to d = dso, assuming that its variation with depth is small. Then, the carbon- 
ion attenuation ratios were determined by reduction of the carbon-ion dose contribution, 
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TABLE II. Theoretical double-ratio factors for water equivalency of plastics. 



Material 


rie / P 

^eui Pw 


S 1 P 


T / S 


a 1 S 


HDPE 


1.027 


1.063 


0.710 


1.110 


PMMA 


0.972 


0.979 


0.890 


1.044 


PET 


0.937 


0.939 


0.953 


1.011 


POM 


0.960 


0.963 


0.969 


1.004 



Ms) = 2M^i2i^, (15) 

which were compared with the linear-attenuation-model calculations. 



no III. RESULTS 



A. Theoretical effective densities 

Table HIl shows the resultant effective-density calculations in forms insensitive to p vari- 
ation. For HDPE of p = 0.96 g/cm^, the electron-density and stopping-power ratios are 
^e/^ew = 0.985 and S/S^, = 1.020, which are close to 1. However, HDPE differs from water 
115 by —29.0% in multiple scattering and +11.0% in nuclear interactions for an equal range 
shift. Similarly, PMMA, PET, and POM differ from water by -11.0%, -4.7%, and -3.1% 
in multiple scattering and by +4.4%, +1.1%, and +0.4% in nuclear interactions. These plas- 
tics are ordered as POM, PET, PMMA, and HDPE by the uniformity of effective densities 
for the relevant interactions, which is highly correlated with the oxygen content. 



120 B. Experimental dose variation 

Figure [2] shows explicit cases for variation of the moderated Bragg peak by insertion of 
range-shifting material, with which the peak and fragment doses were measured. The resul- 
tant measurements are shown in Table IIIIl where D^p and Dq j are the peak and fragment 
doses for s = and available beam range Rq was a sum of the effective thickness of the 
125 tank wall, that of the PPIC cover, and the range-equivalent depth measured for s = 0. For 
proton beams with large range straggling, the 80%-dose depth d^Q is commonly used as the 
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FIG. 2. Variation of moderated Bragg peak by insertion of (a) 8-cm HDPE for 290 MeV and (b) 
20-cm HDPE for 430 MeV (x) with respect to that without an insert (o). 
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FIG. 3. Measured carbon-ion attenuation for range shift by HDPE (x), PMMA (+), and POM 
(o) with hnear-attenuation-model calculations (solid lines). 

range-equivalent depth/^ which may be apphcable to the moderated Bragg peak. By sub- 
tracting the fragment contribution, the range-equivalent depth was redefined to 82%-dose 
depth for Df/Dp = 0.115 and to 84%-dose depth ^§4 for Df/Dp = 0.199. Figure [3] shows 
130 the measurements of a{s) and the linear- attenuation- model calculations of Eq. f|T^ with 
= (0.084 cm-i)/ [(^/a^)/(5/5^)]pMMA = 0-080 cm'^ 




IV. DISCUSSION 

As the density of polymer plastics depends on manufacturing condition, regular calibra- 
tion of the effective density may always be necessary for conversion of thickness to range 
135 shift. Therefore, the stopping-power proximity to water {S ~ Sy^) for HDPE is not an es- 
sential advantage for range control. However, the proximity is desirable when geometrical 
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TABLE III. Results of carbon-ion-beam experiment. 



E/A Ro 

iVlc V Clll 


Material 


t 
cm 


s 

cm 


dso Dp Df 

— a 

cm Dop Dop 


290 14.36 







12.06 1 0.115 1 




HDPE 


4 
8 


4.02 
8.00 


8.04 0.987 0.118 0.982 
4.06 0.974 0.117 0.968 




PMMA 


3 


3.48 


8.58 0.998 0.114 0.999 






7 


8.14 


3.92 0.992 0.113 0.993 




POM 


3 


4.09 


7.97 1.000 0.114 1.001 






6 


8.17 


3.89 1.001 0.114 1.002 


430 28.23 







25.95 1 0.199 1 




HDPE 


10 


9.97 15.98 0.966 0.200 0.956 



20 20.04 5.91 0.937 0.205 0.914 

PMMA 9 10.45 15.50 0.992 0.199 0.990 

17 19.75 6.20 0.981 0.200 0.975 

POM 8 10.88 15.07 1.001 0.198 1.002 

15 20.42 5.53 0.999 0.195 1.003 



properties are of principal interest in dosimetry, such as field width and penumbra size, to 
preserve geometrical effects in water. For applications that allow or require regular air gaps 
in a phantom, higher S/S^ values are acceptable or even preferable. For example, Yajima et 
140 al. developed a multilayer ionization chamber composed of 3-mm and 4-mm PMMA plates 
interleaved with 1-mm air gaps,- which would reduce the density by factor 7/9 resulting in 
mean stopping-power ratio {3/3^) = 0.90. If POM is used instead, the {3/3^) value will 
be 1.06, which will further be adjusted to 1.00 by extending the air gaps to 1.26 mm. 

Beam blurring due to multiple scattering in a range compensator should ideally be min- 
145 imized. In this regard, HDPE with a small {T/T^)/{3/ 3^) value (0.710) is superior to 
PMMA (0.890), PET (0.953), POM (0.969), and water (1), although the effect for carbon 
ions will be small. For example, the rms angle of carbon ions of 20-cm range scattered by 
10-cm water is estimated to be y/{A6'^) = 7.3 mrad by Eq. ([5]). The beam blurring will 
be 1.9 mm for HDPE, 2.1 mm for PMMA, 2.2 mm for PET, and 2.2 mm for POM at the 
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150 isocenter when the range compensator is at 30 cm upstream. 

For nuclear interactions, attenuation of carbon ions in HDPE is excessive. For example, 
the reduction of the Bragg peak was 8.9% with HDPE for 20-cm range shift to a carbon-ion 
beam of 28.23 cm range, while it was 2.5% with PMMA and 0.0% with POM. Invariance 
of the Bragg peak is essential for treatment planning and delivery systems based on water- 

155 phantom dosimetry. A simplistic model estimation based on the Sihver cross-section formula 
and empirical linear- attenuation approximation agreed with the measurement at the level 
of 1%, with which PET is expected to be in between POM and PMMA. 

Material clarity is needed for some dosimetric applications such as water tank. The 
clearest plastic PMMA has been the standard for such purposes. This study suggests that 

160 another clear plastic PET may desirably substitute for PMMA to improve dosimetric water 
equivalency for ion beams. 

V. CONCLUSIONS 

In this study, we evaluated dosimetric water equivalency of four common plastics, HDPE, 
PMMA, PET, and POM, by uniformity of effective densities for carbon-ion-beam interac- 

165 tions. We found that POM was virtually indistinguishable from high-density water, which 
is ideal for range control and preferable for dosimetry with regular air gaps, such as for mul- 
tilayer ionization chamber. For applications that require clarity, PMMA was verified to be 
reasonably water equivalent and PET would be even better. The poorest was HDPE, which 
increased attenuation of carbon ions by 8.9% for 20-cm range shift. HDPE was superior 

170 for small scattering power although the effect would be marginal in practice. We verified 
that water equivalency of plastics was highly correlated with oxygen content. A simplistic 
nuclear-interaction model based on the Sihver cross-section formula and an empirical linear- 
attenuation approximation agreed with the measured attenuation of carbon ions at the level 
of 1%. 
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